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ABSTRACT 

The symbiotic-like binary Mira and nova V407 Cyg was observed in outburst on March 
2010 and monitored in several wavelength bands. The outburst had, to some extent, 
characteristics similar to those observed during other nova eruptions, such as recently 
occurred in RS Oph and U Sco, suggesting that the blast wave interacted with the 
giant companion and propagated through a dense circumstellar medium enveloping 
the binary system. Here we report on multi-dimensional hydrodynamic simulations 
describing the 2010 outburst of V407 Cyg, exploring the first 60 days of evolution. 
The model takes into account thermal conduction (including the effects of heat flux 
saturation) and radiative cooling; the pre-explosion system conditions included the 
companion star and a circumbinary density enhancement that are believed to influence 
the evolution and morphology of the blast wave. The simulations showed that the blast 
and the ejecta distribution are both aspherical due to the inhomogeneous circumstellar 
medium in which they expand; in particular they are significantly collimated in polar 
directions (producing a bipolar shock morphology) if the circumstellar envelope is 
characterized by an equatorial density enhancement. The blast is partially shielded 
by the Mira companion, producing a wake with dense and hot post-shock plasma on 
the rear side of the companion star; most of the X-ray emission produced during the 
evolution of the blast arises from this plasma structure. The observed X-ray lightcurve 
can be reproduced, assuming values of outburst energy and ejected mass similar to 
those of RS Oph and U Sco, if a circumbinary gas density enhancement is included in 
the model. In particular, our "best-fit" model predicts that the 2010 blast propagated 
through a circumbinary gas density enhancement with radius of the order of 40 AU 
and gas density ss 10 6 cm' 3 and that the mass of ejecta in the outburst was M e j « 
2 x 10~ 7 Mq with an explosion energy Eq ps 2 x 10 44 erg. Alternatively, the model 
can produce a similar X-ray lightcurve without the need of a circumbinary gas density 
enhancement only if the outburst energy and ejected mass were similar to those at 
the upper end of ranges for classical novae, namely M j w 5 x 10~ 5 M Q and E « 
5 x 10 46 erg. 

Key words: shock waves - binaries: symbiotic - circumstellar matter - stars: indi- 
vidual (V407 Cyg) - novae, cataclysmic variables - X-rays: binaries 



1 INTRODUCTION 

The symbiotic binar y V407 Cyg was discove red in outburst 
at 7th magnitude bv lNishivama et al.l |2010l ) on 2010 March 
10, in stark contrast to occasional brightenings observed 
over the years by a few magnitudes compared to histor- 
ical low states in the range 13-17m pg . The binary com- 
prises a white dwarf (WD) and Mira-type M6-7 red gi- 



ant pulsating with a period of 76 3 d, with an orbital pe - 
riod of 43 yr dMunari et al.l Il99d: iKolotiloy et alj l2003l h 
Spectroscopy obtained bv lMunari et al.l |201ll 'l and C. BuiQ 
two days after outb urst showed very broad Balmer lines; 
iMunari et alj l^Oll! - ) reported a full-width at half-maximum 
(FWHM) of 2760 km s" 1 on day 2.3, and described the 
outburst as a He/N nova expanding within the wind of the 
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Mir a companion, similar to the 2006 explosion of RS Oph 



(see A. Evans, M. F. Bode. T. J. O'Brien, fc M. J. Darnlevl 



I200E . and references therein). SiO maser emission was 
found to dec l ine d ramatically after the burst, leading 
iDeguchi et alj (1201 ll ) to conclude that the maser emitting 
regions were essentially wiped out in a time scale of two 
weeks by the propagation of the nova shock into the Mira 
envelope. 

The RS Oph blast wave produced copious X-rays from 
at most 3 days after optical discovery, when it was first 
observed by Swift and RXTE, characterized by a steady 
decline thereafter dSokoloski et all IgOOfj ; iBode et all 120061 ; 
iNelson et al.ll20Qg| ; IPrake et al.ll2009l ). V407 Cyg also devel- 
oped into an X-ray source, though somewhat more slowly, 
exhibiting a strong brightening at 20 days o r so, a peak near 
30 days and a steady decline. I Shore et alj (2011) analyzed 
the Swift UV and X-ray lightcurves together with optical 
spectra obtained during the first three months of the erup- 
tion. B aimer lines s howed steady secular narrowing (see also 
iMunari et al]|201ll ). and lines of [Ca V], [Fe VII], [Fe X], and 
He II exhibited asymmetric profiles attributed to an aspher- 
ical expansion remi niscent of that dia gnosed for RS Oph 
from X-ray lines bv lDrake et~al] (120091 ). 

While not nearly as bright as RS Oph in X-rays, 
V407 Cyg was in some ways even more remarkable, with a 
firm detection by the Fermi Large Area Telescope of variable 
7-ray emission in the 0.1-10 GeV range on 2010 March 10 — 
the sa me day as optica l discovery — that persisted for two 
weeks l|Abdo et al.ll20ld) . This was the first ever 7-ray detec- 
tion o f a nova explosion and was attributed by lAbd o et alj 
<|2010l ) to 7T° production and subsequent decay resulting 
from c ollisions of protons accelerated in the shock. iLii et all 
(2011) pointed out that short period symbiotic novae such as 
RS Oph are unlikely to produce 7-rays because of the fast 
evolution of the blast wave in the higher density environ- 
ment closer in to the red giant, whereas long period systems 
such as V407 Cyg pro vide a much longer-lived accelerator. 
iRazzaaue et alj (|2010h further noted that V407 Cyg could 
have been a detectable source of high energy neutrinos. 

The outburst mechanism for novae is thermonuclear 
runaway on the WD triggered by the mass of accreted 
mate rial exceeding a critical limit (|Starrfield et alj 1 19851 . 
1 19881 ). In two earlier papers we have applied sophisticated 
hydrodynamic mode ls to the nov a explosions on RS Oph 
<|Orlando et alj 120091 ) and U Sco (|Drake fc Orlandol |2010| ) 
to provide key insights into the nature of the explosion and 
its environment. In both sets of models, circumbinary mate- 
rial proved crucial in explaining the observed X-ray emission 
while providing a degree of collimation to the explosions. 
Here we describe similar detailed hydrodynamic simulations 
of the V407 Cyg explosion, and pay particular attention to 
the effects of circumbinary gas density enhancement and the 
secondary Mira companion on the explosion. In Sect. [2] we 
describe the hydrodynamic model, the numerical setup, and 
the synthesis of X-ray emission; in Sect. we discuss the 
results; and finally in Sect. [4] we draw our conclusions. 



2 HYDRODYNAMIC MODELING 

The blast wave was modeled by numerically solving the 
time-dependent fluid equations of mass, momentum, and en- 



ergy conservation, including radiative losses described by an 
optically-thin plasma and thermal conduction; the latter in- 
corporated the effects of heat flux saturation. Owing to the 
long timescale of the blast evolution, radiative losses and 
thermal conduction were found to be more im portant than 
in the cases of RS Oph (|Orlando et all [2009) and U Sco 
l|Drake fc Orlandol [2010) . The long evolution timescale of 
the blast coupled with large expansion velocities of a few 
thousand km s~* rendered the spatial extent too large and 
computationally demanding to perform an extensive set of 
fully 3-dimensional (3D) hydrodynamic simulations while 
still resolving the structure of the immediate binary envi- 
ronment, even with many levels of adaptive mesh refine- 
ment. The computational cost is raised by the inclusion of 
thermal conduction that is solved explicitly, such a scheme 
being subject to a rather restrictive stability condition as the 
thermal conduction timescale is generally shorter than the 
dynamical one (e.g. jHujeirat fc Camenzind 2000; H uieirad 
120051 ; lOrlando et al.ll2f7o8l , l2010l ). 

Given the large computational cost required by 3D sim- 
ulations, we adopted the following strategy: we first explored 
the wide parameter space of the model by adopting a 2- 
dimensional (2D) cylindrical coordinate system (r, z) and, 
therefore, assuming the system to be symmetrical with re- 
spect to the axis passing through the WD and the com- 
panion star; then, for the set of parameters found to best 
reproduce the observations, we relaxed the hypothesis of 
axisymmetry and performed a fully 3D simulation in carte- 
sian geometry. The 2D and 3D mode ls developed h e re are 
othe rwise similar to th e 3D m odels of lOrlando et al.l (|2009l ) 
and IPrake fc Orlandd (|2010l ) and we refer to those works 
for further details. The calculations were performed using 
FLASH, an adaptive m esh refinement mult iphysics code for 
astrophysical plasmas (|Frvxell et al 1 120001 ) extended with 
additional computational m odules to handle radi ative losses 
and thermal conduction (see lOrlando et aT1 l2005 for the de- 
tails of the implementation). The hydrodynamic equations 
for compressible gas dynamics are solved using the FLASH 
implementation of the pie cewice-parabolic method (PPM; 
IColella fc Woodward! Il983 ). 

W e adopted the system parameters of IMunari et al.l 
(1990); these are listed in Table [1] For the hydrodynamic 
models, the most important system parameters are the or- 
bit semi-major axis, the wind mass loss rate and wind ter- 
minal velocity. These determine the gas density into which 
the blast occurs, and for a given explosion energy largely 
control the subsequent evolution timescale of the resulting 
shock wave system. The mass lo ss rate and term i nal ve - 
locity are based on the data of iKnapp fc Morris! (|l985h . 
while the orbital separation was estimated based on dust 
extinc tion changes assume d connected with orbital modu- 
lation. IMunari et al.l (|l990l ) estimated orbital separations of 
14.0, 15.5 and 16.4 AU for WD masses M W d = 0.5, 1.0 
and 1.4 Mq, respectively; we adopt 15.5 AU and Mwd = 
1.0 Mq. We also adopted t he IMunari et all (1990) distance 
of 2.7 kpc, derived from the lGlass fc FeastJ (|l982l T luminosity 
vs. pulsation period relation extrapolated to the V407 Cyg 
period of 745 days. 

Figure [T] shows an example of the initial conditions 
adopted in our simulations. The thermonuclea r explosion i s 
initiated by a spherical Sedov-type blast wave (|Sedovlll959l ) 
centered on the WD, with radius ro = 10 8 km (red circle 
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Table 1. Adopted parameters and initial conditions for the hydrodynamic models of the 2010 V407 Cyg explosion 
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Parameter 



Value 



Secondary star radius 
Binary separation 
Orbital period 
Spatial domain 



AMR max. resolution 
Time covered 



R cs = 2.2 AU 
a = 15.5 AU 
P = 43 yr 

(2D cylindrical) (3D cartesian) 

< r < 530 AU -530 < x < 530 AU 

-530 s£ z s£ 530 AU ^ y ^ 530 AU 

0^ 530 AU 
5 x 10 11 cm (3.33 x 10 -2 AU) 
0-60 days 



Model abbreviation 


geometry 


Eo 
[erg] 


[Mq] 




[cm -3 ] 


^cde 

[cm -3 ] 


h 

[AU] 


h 

[AU] 


h 

[AU] 


E44.3-NW7 


2D 


cylindrical 


2 X 10 44 


2 x 10" 


7 


10 7 










E44.3-NW8 


2D 


cylindrical 


2 X 10 44 


2 x 10- 


7 


10 8 










E44.3-NW9 


2D 


cylindrical 


2 x 10 44 


2 x 10" 


7 


10" 










E45-NW8 


2D 


cylindrical 


10 45 


10~ 6 




10 s 










E45-NW9 


2D 


cylindrical 


10 45 


10" 6 




10 9 










E46-NW9 


2D 


cylindrical 


10 46 


10" 5 




10 9 










E46.6-NW9 


2D 


cylindrical 


5 x 10 46 


5 x 10" 


5 


10 9 










E47-NW10 


2D 


cylindrical 


10 47 


10- 4 




10 10 










E44-NW7-CDE6-L40 


2D 


cylindrical 


10 44 


io- 7 




10 7 


10 6 


40 


40 




E44-NW7-CDE6.7-L40 


2D 


cylindrical 


10 44 


lO" 7 




10 7 


5 x 10 6 


40 


40 




E44-NW7-CDE7-L40 


2D 


cylindrical 


10 44 


io- 7 




10 7 


10 7 


40 


40 




E44.3-NW7-CDE6.3-L20 


2D 


cylindrical 


2 x 10 44 


2 x 10- 


7 


10 7 


2 x 10 6 


20 


20 




E44.3-NW7-CDE6.3-L40 


2D 


cylindrical 


2 x 10 44 


2 x 10- 


7 


10 7 


2 x 10 6 


40 


40 




E44.3-NW7-CDE6.3-L80 


2D 


cylindrical 


2 x 10 44 


2 x 10" 


7 


10 7 


2 x 10 6 


80 


80 




E44.7-NW7-CDE6.7-L40 


2D 


cylindrical 


5 x 10 44 


5 x 10" 


7 


10 7 


5 x 10 6 


40 


40 




E45-NW7-CDE6-L40 


2D 


cylindrical 


10 45 


10~ 6 




10 7 


10 6 


40 


40 




E45-NW7-CDE6.7-L40 


2D 


cylindrical 


10 45 


IO" 6 




10 7 


5 x 10 6 


40 


40 




E45-NW7-CDE7-L40 


2D 


cylindrical 


10 45 


IO" 6 




10 7 


10 7 


40 


40 




3D-E44.3-NW7-CDE6.3 


3D 


cartesian 


2 x 10 44 


2 x 10- 


7 


10 7 


2 x 10 6 


53 


53 


27 



on the right in the inset panel of Fig. HJ. In analogy with 
supernova explosions, the total energy of the blast, Eo, is 
partitioned so that most of the explosion energy is kinetic 
(in particular we assumed 1/4 of the energy contained in 
thermal energ y and the other 3/4 in kinetic energy as repre- 
sentative; e.g. iDohm-Palmer fc Joneslll996h . The total mass 
of the ejecta is M e j. The blast propagates through the ex- 
tended outer atmosphere (the wind) of the Mira companion 
and is off-set from the origin of the wind density distribution 
by 15.5 AU (i.e. the system orbital separation; see inset panel 
of Fig. [TJ. We assumed the gas density in the wind is propor- 
tional to R~ 2 (where R is the radial distance from the Mira). 
In addition to the R~ 2 density distribution, we also included 
a circumstellar envelope (or a circumbinary density enhance- 
ment in the wind, hereafter CPE), as predicted by detailed 
hydrodynamic mod eling (e.g. iMastrodemos fc Morris! 1 19991 ; 
Walder et al.| [2008). In particular, in systems comprising a 
red giant star, these models predict a CDE created by grav- 
itational accumulation of the cool red giant wind toward the 
WD and a spiral shock wave caused by the motion of the 
stars through this gas. This scenario is also supported by 
observations in different wavelength b ands of similar nova 
outbursts as, for i nstance, RS Oph (e .g. O'Brien et al. 2006J ; 
iBode et alj|2007l : iRupen et all [20(53 ; brake et al.ll2009l ). In 



P = p™ 



+ Pcde e 



[-(x/h) 2 -{y/h) 2 -(z/h) 2 ] 



(1) 



cartesian geometry, the mass density distribution of the un- 
perturbed circumstellar medium (CSM) is given by: 



where p w = [imnriv, is the wind mass density close to the 
surface of the Mira companion, p = 1.3 is the mean atomic 
mass (assuming cosmic abundances), bih is the mass of the 
hydrogen atom, R ca is the radius of the Mira, R is the ra- 
dial distance from the Mira, p c d c = /iWHicdc is the density 
enhancement close to the Mira, and Zi, h and I3 are char- 
acteristic length scales determining the size and shape of 
the CDE. As examples, Fig. [2] shows the density profiles 
along the x axis for a model including the CDE (run 3D- 
E44.3-NW7-CDE6.3) and for a model without a CDE (run 
E44.3-NW7). The Mira companion is included in the calcu- 
lation as an impenetrable body with radius R cs = 2.2 AU 
(i.e. « 500 Rq; white circle on the left in the inset panel of 
Fig. El. 

In the 3D simulation, the hydrodynamic equations were 
solved in one quadrant of the whole spatial domain with the 
coordinate system oriented in such a way that both the WD 
and the companion star lie on the x axis (see Fig. [1] see also 
iDrake fc Orlandoll2010l for more details). The companion is 
at the origin of the coordinate system, (x,y,z) = (0,0,0), 
and the computational domain extends 1060 AU in the x 
direction and 530 AU in both the y and z directions; the 
WD is located to the right on the x axis (y = z = 0) at 
x — 15.5 AU. We imposed reflecting boundary conditions at 
J/min = and Zmi n = (consistent with the adopted sym- 
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metry) and outflow (zero-gradient) conditions at the other 
boundaries. 

Our 2D simulations considered a slab corresponding to 
the (x, z) plane of the 3D simulation, and the hydrodynamic 
equations were solved in one half of the 2D spatial domain, 
taking advantage of the symmetry of the configuration. In 
particular, assuming axisymmetry about the axis passing 
through the companion star and the WD, the cylindrical 
coordinate system was oriented in such a way that both the 
donor star and the WD lie on the z axis, with the former at 
the origin of the coordinate system, (r,z) = (0,0), and the 
latter at z = 15.5 AU. Since both the donor star and the 
initial nova remnant lie on the symmetry axis they are mod- 
elled as spheres in the cylindrical coordinate system. The 
computational domain extended 530 AU in the r direction, 
and 1060 AU in the z direction. The adopted geometrical 
system configuration implied imposition of an axisymmet- 
ric (reflecting) boundary condition at r m i n = and outflow 
conditions at the remaining boundaries. 

In order to compare the results of the 2D simulations 
with those derived from the 3D simulation, we reconstructed 
the 3D spatial distributions of all the physical variables (e.g. 
mass density and temperature) derived with the 2D simu- 
lations, by rotating the 2D slab about the z axis, according 
to the symmetry of the problem. Thus we passed from the 
cylindrical coordinate system to the cartesian one, orienting 
the cartesian coordinate system in such a way that both the 
WD and the companion star lie on the x axis, as in the 3D 
simulation. 

The explosion and subsequent blast wave was followed 
for a total of 60 days in order to explore the evolution of the 
X-ray emission and study the effects of the circumstellar en- 
vironment on the evolution of the blast during the evolution- 
ary phase that was characterized by a pronounced pea k of 
X-ray emission in the 2010 outburst (jShore et al.ll201lT ). As 
with the modeling of the U Sco blast, the small scale of the 
stellar system compared with the size of the rapidly expand- 
ing blast wave over the 60 day period of interest presents a 
major computational challenge. To this end, we exploited 
the adaptive mesh capabilities of the FLASH code by using 
12 nested levels of adaptive mesh refinement, with resolution 
increasing twice at each refinement level. This grid configu- 
ration yielded an effective maximum resolution of « 5 x 10 6 
km at the finest level, corresponding to « 20 grid points per 
initial radius of the blast and » 66 grid points per radius of 
the Mira companion. The calculations were performed using 
an automatic mesh derefinement scheme in the whole spatial 
domain except in the Mira (where the resolution does not 
vary during the blast evolution) that kept the computational 
cost approximately constant as the blast expanded. The ef- 
fective mesh size was 16384 x 32768 for the 2D simulations, 
and 32768 x 16384 x 1638 4 for the 3D simulat i on. 

Solar abundances of iGrevesse fe Sauva 3 |l998t ) (GS) 
were assumed for the wind and circumstellar density en- 
hancement, while ejecta metal abundances were assumed 
to be enhanced by a factor of t en. This latt e r cho ice was 
guided by the Swift spectra that IShore et al.l l |201ll ) found 
to be rich in both N and O — by possibly more than a fac- 
tor of ten — and by ob servations of He-rich ej ecta in the 
out burst of U S c o (see iDrake fc Orlando! l2010h . as well as 
the IDrake et all <|20Qgf > high-resolution X-ray spectroscopic 
study of the 2006 RS Oph blast that found evidence for 
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Figure 1. Colour-coded cross-section image of the gas density 
distribution, on a logarithmic scale, in units of cm -3 , showing the 
initial conditions of run 3D-E44.3-NW7-CDE6.3. The inset panel 
shows a close-up view of the initial geometry of the V407 Cyg 
system. The companion star is at the origin (white circle on the 
left of the inset), and the white dwarf lies on the x axis at x = 15.5 
AU (red circle on the right). 

metal-rich ejecta. The adopted abundances are relevant for 
the radiative losses from shocked ejecta, and for the lo- 
cal absorption by the shocked CSM (with GS abundances) 
and by the ejecta (with GS abundances x 10) encountered 
within the blast wave. The choice of abundances is also rel- 
evant for computation of the emitted X-ray intensity of the 
blast. Such emission was synthesized from the mode l result s 
using the methodology described bv lOrlando et all {2009). 
The synthesis includes: thermal broadening of emission lines, 
Doppler shift of lines due to the component of plasma veloc- 
ity along the line-of-sight, photoelectric absorption by the 
interstellar medium (ISM), CSM, and ejecta. The absorp- 
tion by the ISM was calculated assuming a column density 
Ah = 2 x 10 21 cm -2 ; the local absorption was calculated 
self-consistently from the distributions of shocked CSM and 
shocked ejecta. 

The influence of the different system parameters was in- 
vestigated through the 2D simulations, by exploring models 
with an initial energy of explosion, Eo, in the range 10 44 - 
10 47 erg, ejecta mass, M e j, in the range 10 _7 -10 -4 Mq, 
wind density in the range 10 7 -10 10 cm -3 , CDE density in 
the range 10 6 -10 7 cm 3 (see Table [TJ. The ranges of ejected 
mass and outburst energy include those typical of recurrent 
novae (10" 7 < Af cj < 10~ 6 Mq, 10 44 < E < 10 45 erg), 
given the close similarities between the 2010 outburst of 
V407 Cyg and those of RS Oph and U Sco, and the ranges 
typical of classical novae (10~ 5 < M ej < 10" 4 Mq, 10 46 < 
E < 10 47 erg; e.g. iBoddl^Olol '). In the 2D simulations, the 
CDE has a spherical shape (according to the symmetry of 
the system) and Eq. [T] reduces to 

P-P-(^) 2 +P^ el ~ (r/h)2 ' (Z/l2)2] ( 2 ) 
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Figure 2. Initial density profiles along the x axis for a model 
including the CDE (run 3D-E44.3-NW7-CDE6.3; solid line) and 
for a model without CDE (run E44.3-NW7; dashed line). The 
dotted lines mark the radius of the Mira companion, R cs , and 
the binary separation, a. 



with l\ = h- The CDE thickness was explored in the range 
20-80 AU. The additional 3D simulation allowed us to study 
the collimation of the ejecta due to an equatorial density 
enhancement: in this case we adopted the parameters of the 
2D simulation that best fit the observati ons and explored the 
effects of a disk-like CDE (according to lWalder et alj|2008h 
on the blast evolution, by assuming fa = h and Z3 < fa in 
Eq.DG 



3 RESULTS AND DISCUSSION 
3.1 Hydrodynamic evolution 

The V407 Cyg models show a blast wave evolution that is in 
so me ways intermediate between those pre sented for RS Oph 
bvlOrlando et all (|2009l ) and for U Sco bv lDrake fc Orlandol 
(2010). The WD of V407 Cyg is characterized by the massive 
circumstellar gas envelope of its Mira companion in which 
it sits, similar to the WD of RS Oph located within the 
dense wind of its red giant secondary, with no evident well- 
defined accretion disk (unlike U Sco) . V407 Cyg has a much 
wider orbital separation than RS Oph and U Sco, implying 
that the nova explodes initially into a relatively low density 
environment, as in the case of U Sco, despite the presence of 
the CDE. The size of the Mira companion (R cs = 2.2 AU) 
is significant compared to the size of the system, so that the 
Mira partially shields the blast, as found in U Sco. 

The evolution of the blast for three representative 
V407 Cyg models is illustrated in Fig. [3] Post-blast gas den- 
sity distributions in the (x, z) plane bisecting the system 
(the plane of the orbital axis; the equatorial plane is edge 
on) are shown, sampled at the labeled times. The figure 
also shows the distribution of ejecta within the blast wave 
by highlighting regions where more than 90% of the mass 
is material ejected in the explosion (dashed contours). The 
hot post-shock plasma is dominated by thermal conduction 
rather than radiative cooling, and this partially suppresses 
the hydrodynamic instability that would oth erwise develop 
during the evolution of the blast wave (see lOrlando et al.l 



I2005I .I 2008). Owing to the low density environment in which 
the blast propagates, and to the high post-shock tempera- 
ture, there are in fact no regions within the blast that appear 
dominated by radiative losses. Thermal conduction is ex- 
pected to dominate over radiative cooling when the conduc- 
tion time-scale is shorter than the cooling time-scale or, in 
other words, in plasma structures with a length-scale smaller 
than the Field length, Lf ~ 10 6 T 2 /n H (|Begelman fc McKeei 
Il99d ; lOrlando et ai1l2005l ), where T and hh are the plasma 
temperature and particle number density, respectively. For 
our models, Lp ~ 240 AU, such that this length-scale is 
larger than the size of the plasma structures produced in 
the blast. 

The models all show aspherical shock morphologies, 
rendered by the blast wave propagating through the inho- 
mogeneous circumstellar gas distribution. In models without 
a CDE (hereafter NOCDE models; e.g. run E44.3-NW8 in 
Fig. [3}, the aspherical morphology is mainly due to the off- 
set of the blast wave with respect to the origin of the wind 
density distribution. A similar morphology is found in mod- 
els including a spherical CDE (hereafter SPHCDE models) 
with an origin coincident with that of the wind density dis- 
tribution (e.g. run E44.3-NW7-CDE6.3-L40 in Fig. [3$. Here 
again the morphology is determined by the off-set of the 
blast wave with respect to the origin of the CDE. At vari- 
ance with NOCDE models, however, the blast initially prop- 
agates in a medium with a density distribution decreasing 
away from the companion star more slowly than the R~ 2 
wind density distribution (see Fig. 

The CDE has quite an important effect on the blast 
evolution if its shape is disk-like, as in run 3D-E44.3-NW7- 
CDE6.3 (see right panels in Fig- [3jl . In this case, the blast 
morphology is similar t o that found f rom the modeling of 
the R S Oph outburst (|Walder et all 120081 ; lOrlando et all 
2009): both the blast and ejecta are strongly collimated in 
polar directions, leading to a bipolar shock morphology. It 
is remarkable that the shock expansion at day 30 is almost 
parallel to the z-axis. Figure [4] shows, as an example, the 
collimation of blast and ejecta 60 days after the outburst 
for the model 3D-E44.3-NW7-CDE6.3. A similar poleward 
blast collimation was a feat ure of hydrodynamic si mulations 
of the early U Sco blast (|Drake fc Orlandol [2oToh : there, 
the collimation was caused by a dense accretion disk sur- 
rounding the WD. Such a blast collimation is also suggested 
by observations. Indeed, collimation signatures found during 
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represent a growing consensus in the literature that blast 
collimation is a common feature of nova outbursts. 

In all the models, the shock front propagating toward 
the Mira companion is partially shielded by it and refracted 
around it. As a result, the shock converges on the rear side 
of the companion star, undergoing a conical self-reflection 
there, and producing a wake with relatively dense (nn ~ 
10 6 cm" 3 ) and hot (T » 5 x 10 7 K) post-shock plasma. We 
expect, therefore, that this structure may lead to conspic- 
uous observable X-ray emission. The potential for such a 
structure to explain persistent blue- war d extension of opt i- 
cal emission lines was also suggested bv lShore et al.l (|201lh . 
At early stages, the secondary star also gives rise to a bow 
shock with a temperature T « 10 7 K. However, this feature 
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Figure 3. Cross-section images of the gas density distribution, on a logarithmic scale in units of cm -3 (see the color table in the upper 
right corner of each panel), sampled at the labeled times for three representative V407 Cyg models. The Mira companion is at the origin 
and the WD is off-set from the origin to the right by 15.5 AU. The different panels illustrate the evolution of the blast for a model 
without CDE (run E44.3-NW8; left panels), for a model with a spherical CDE (run E44.3-NW7-CDE6.3-L40; center panels), and for a 
model with a disk-like CDE (run 3D-E44.3-NW7-CDE6.3; right panels). Inset panels show the blast structure closer to the system origin. 
The white circle in the inset represents the Mira companion. The white dashed contour encloses the ejecta. The white solid contour 
denotes the regions with plasma temperature T > 10 6 K. 



is characterized by relatively low densities (jih » 10 5 cm -3 ; 
one order of magnitude lower than the density in the wake) 
and we expect only a modest contribution to X-ray emission. 



3.2 X-ray emission 

From the simulations, we synthesized the X-ray emission in 
the [0.6 — 12.4 ] keV band, using the method described by 
lOrlando et ail (|2009D and outlined in Sect. [2] Fig. [5] shows 
the maps of X-ray emission projected along the line of sight 



at day 30 (namely when a pronounced peak of X-ray emis- 
sion was observed in the 2010 outburst) for the models re- 
ported in Fig. [3] The maps are shown on linear (upper pan- 
els) and logarithmic (lower panels) scales to highlight struc- 
tures with very different emission levels. The plane of the 
orbit of the central binary system lies on the (x, y) plane 
and is assumed to be edge on. In all cases, most of the X-ray 
emission arises from the high-temperature shocked CSM in 
the wake produced by the convergent shock on the rear side 
of the companion star. Similar plasma structures have been 
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Figure 4. Three-dimensional rendering of particle number den- 
sity, on a logarithmic scale, in units of cm -3 (see color table in 
the upper left corner of the panel), at day 60 for the model 3D- 
E44.3-NW7-CDE6.3. The ejecta are highlighted in red. The plane 
of the orbit of the central binary system lies on the (x,y) plane. 



found in nu merical simulations describing the outburst of 
RS Oph (e.g. lWalder et al . 2008) and are generally expected 
as a result of the shielding of the blast by the companion 
star. 

The X-ray lightcurves of the simulations in the [0.6 — 
12.4] keV band are shown in Fig. [6] together with the ob- 
served lightcurve. The l atter was based on the <SVre/£/XRT 
count rates reported by IShore et al.l (|201ll ). considering a 
maximum X-ray luminosi ty of Lx = 1.2 x 10 34 erg s" 1 
at a dist ance of 2.7 kpc JShore et al.l [20111 ). As found for 
U Sco bv lDrake fc Qrlandol l|2010l ). in NOCDE models (left 
panel in Fig. [6]) in general the X-ray luminosity, Lx, is 
larger for higher explosion energy (and ejecta mass), and 
wind gas density. For outburst energies and ejecta masses 
in the range of values typically observed for recurrent novae 
(i.e. lO" 7 < M ej < 10" 6 M Q , 10 44 < E < 10 45 erg; e.g. 
lBodell2010l ), the synthetic Lx reaches its maximum very 
quickly, within few days after the explosion, so that none 
of these lightcurves fit the observed one which is charac- 
terized by a maximum around day 30. On the other hand, 
assuming outburst energies and ejecta masses typical of clas- 
sical novae (i.e. in the ran ges 10~ 5 < M c j < 10~ 4 Mq, 



10 46 < E < 10 47 erg; e.g. iBodd lioiO). the peak of X-ray 



emission can be reached at later times. Among these mod- 
els, that which best reproduces the X-ray peak at day 30 is 
E46.6-NW9 (orange line in the left panel in Fig. [6]), although 
it fails in reproducing the descending slope of the observed 
lightcurve. This model predicts outburst energy of the order 
of w 5 x 10 46 erg and mass of ejecta « 5 x 10~ 5 Mq that 
seem to be too high in the case of V407 Cyg that shows 
many similarities with the novae RS Oph and U Sco. 

In the presence of a CDE, the X-ray luminosity mainly 
depends on the physical characteristics of the density en- 
hancement: the greater the size and gas density of the CDE, 
the larger Lx (central and right panels in Fig. [5]). In these 
models, the blast reaches its maximum X-ray luminosity 



later than in NOCDE models, and the time of maximum 
depends on the size and gas density of the CDE and on the 
explosion energy (and ejecta mass). In particular, the time 
delay of maximum X-ray luminosity is longer for lower gas 
density and greater size of the CDE, and for higher explosion 
energy. Among the SPHCDE models, that which best repro- 
duces the observations is E44.3-NW7-CDE6.3-L40 (blue line 
in the right panel of Fig. [6]): it describes a continuous rise of 
X-ray luminosity up to a maximum of Lx ~ 2 x 10 34 erg s _1 
at day 30 and then a subsequent decay phase until day 60, 
with the same slope as the observed lightcurve. It is worth 
noting that this model assumes outburst energies and ejecta 
masses in the range of values typically observed for recurrent 
novae, while the "best-fit" NOCDE model assumes values 
typical of classical novae. For E44.3-NW7-CDE6.3-L40, we 
extended the simulation in order to cover the whole period 
in which the observed lightcurve is defined. We found that, 
after day 60, the synthetic lightcurve fades faster than the 
observations which seem to be characterized by a plateau 
between days 60 and 90. This discrepancy could be ex- 
plained as being due to the details of the density struc- 
ture of the CDE; for instance, a radial density profile of 
the CDE decreasing away from the companion star slower 
than that modelled here may change the descending slope 
of the lightcurve. Another possible cause of the discrepancy 
may be some departures from the simple wind density profile 
adopted here. 

The effect of the shape of the CDE on the synthetic X- 
ray image and lightcurve was explored by performing a 3D 
simulation describing a disk-like CDE (run 3D-E44.3-NW7- 
CDE6.3; see right panels in Fig. [3]) with the same parameters 
of run E44.3-NW7-CDE6.3-L40 but with h = h = 53 and 
^3 = 27 in Eq.[T](see Table[T]). In this case, again most of the 
X-ray emission originates in the high-temperature shocked 
plasma in the wake of the companion star (see right panels 
in Fig. 0), although the X-ray source here is more compact 
than in run E44.3-NW7-CDE6.3-L40. The X-ray lightcurve 
derived from run 3D-E44.3-NW7-CDE6.3 is illustrated in 
the right panel of Fig. [(J] (magenta line): the curve is similar 
to that derived from run E44.3-NW7-CDE6.3-L40 and is in 
good agreement with the observed lightcurve. 



4 SUMMARY AND CONCLUSIONS 

We have modeled the first 60 days of evolution of the 2010 
blast wave of the nova V407 Cyg. The models explore the 
3D structure of the blast wave originating from the nova 
outburst, taking into account simultaneously the radiative 
cooling and thermal conduction (including heat flux satura- 
tion). Our analysis indicates the following: 

(i) the blast wave and the ejecta distribution are both 
aspherical due to the inhomogeneous CSM in which the blast 
expands; 

(ii) in case of a disk-like CDE, the blast and the ejecta 
are collimated in polar directions, leading to a bipolar shock 
morphology; 

(iii) the shock front propagating toward the Mira is par- 
tially shielded by it, producing a wake with dense and hot 
post-shock plasma on the rear side of the companion star; 

(iv) during the evolution of the blast, most of the X-ray 
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Figure 5. X-ray images (normalized to the maximum of each panel) in [0.6 — 12.4] keV band projected along the line of sight of the 
blast after 30 days of evolution, corresponding to the density distributions in the middle panels in Fig. \3\ Upper panels are rendered 
with a linear surface brightness scale and lower panels with a logarithmic scale. Both the Mira companion and the WD lie on the x-axis. 
The plane of the orbit is edge on. Most of the X-ray emission originates from the shocked CSM. The white dashed contour encloses the 
ejecta. The white solid contour denotes the forward shock. 




time [ days ] time [ days ] time [ days ] 

Figure 6. X-ray lightcurves in the [0.6 — 12.4] keV band originating from the nova outburst (colour-coded solid lines) for the models 
reported in Table [T] The left panel shows the results for the NOCDE models. The X-ray flux is obtained by integrating the synthetic 
X-ray emission in the whole spatial domain and takes into acco unt the effects of local absorption by shocked CSM and ejecta. The 
reference black dashed line is the observed lightcurve reported bv lShore et al. 
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emission arises from the dense and hot post-shock plasma 
in the wake of the companion star; 

(v) the observed X-ray lightcurve can be reproduced 
without including the CDE in the model only if the out- 
burst energy and ejected mass were of the same order of 
magnitude as those at the upper end of ranges for classical 
novae, namely M e j « 5 x 1CT 5 M Q and _Eo ~ 5 x 10 46 erg, 
i.e. values much higher than those estimated in the case of 
RS Oph and U Sco; 

(vi) the model best reproducing the observations assumes 
the presence of a CDE and suggests that the explosion en- 
ergy in the 2010 outburst was ~ 2 x 10 44 erg (and the mass 
of ejecta ~ 2 x 10~ 7 Mq, i.e. values comparable to those 
estimated for RS Oph and U Sco) and that the blast propa- 
gated through a CDE with radius of the order of 40 AU and 
gas density « 2 x 10 6 cm~ 3 . 

IShore et af] l|201ll ) note that the 2010 outburst of 
V407 Cyg closely resembles the spectroscopic development 
of RS Oph which is considered the prototype of symbiotic- 
like recurrent novae. This fact has suggested that V407 Cyg 
could belong to the same class of recurrent novae, although 
its outburst was the first recorded no va event for this sys tem. 
On the other hand, as discussed by IShore et al.l (|201lT ) the 
historical lightcurve of V407 Cyg is characterized by many 
gaps and it is possible that outbursts similar to that ob- 
served in 2010 were missed. Munari et al (1990) estimated 
an accretion rate of ~ 1O _8 M0 yr _1 for a WD mass of 
Mwd ~ IMq. If this accreted mass were ejected in the 
nova event, our estimate of the ejecta mass for the CDE 
case implies a recurrence time for nova events of 20 years or 
so. For this accretion rate, the recurrence timescale is consis- 
tent witlijii^he^jna^wlrite dwarfs (e.g. Mwd ~ 1.25 M©; 
e.g. iPrialnik fc Kovet j[T995l V Within the fairly large uncer- 
tainties of accretion rate and ejecta mass estimates, this is 
consistent with either a scenario in which one or more out- 
bursts have been missed, or with a recurrence timescale of 
the order of 100 years. For a low mass WD the recurrence 
timescale is two orders of magnitude longer. 

In the light of the above considerations, and the results 
of numerical models showing that gas does accumulate in 
the orbital plane of a symbiotic binar^Q, we suggest that the 
model best describing the 2010 outburst of V407 Cyg is that 
assuming the presence of a CDE and predicting outburst 
energy and mass of ejecta similar to those estimated for 
RS Oph (run 3D-E44.3-NW7- CDE6.3). 

Finally, the X-ray structure predicted by our model has 
interesting implications for the line profiles expected to be 
observed during the outburst. Simulations of the 2006 out- 
burst of RS Oph, that occurred close to orbital quadra- 
ture, have shown that most of the X-ray emission arises 
from a compact source propagating perpendicular to the 
line of sight, which is similar to that found here. This 
source produced asymmetric and slightly blue-shifted line 
profiles due to substantial X-ray absorption of red-shifted 

2 The gravitational accumulation of gas toward the WD is ex- 
pected to have for V4 7 Cyg a form similar to that for RS Oph. 
In fact, IWalder ct al. (2008) have characterized the density en- 
hancement by the ratio R = Vrg / Vorb (where Vrg is the termi- 
nal wind velocity and V or b = 2ira/P, with P the orbital period) 
which for both RS Oph and V407 Cyg is of order 1. 



emission by shocked material |Orlando et alj|2009h . These 
m odel predicti o ns we re in nice agreement with the findings 
of iDrake et al.l (2009) who reported similar asymmetric line 
profiles from the analysis of Chandra/HETG observations of 
RS Oph. In the light of the above results, we expect therefore 
that even in the case of V407 Cyg the line profiles should 
be asymmetric and blue-shifted due to X-ray absorption 
of red-shifted emi ssion by shocked material. Interestingly, 
I Shore et al.l l|201ll ) reported asymmetric line profiles with 
extended blue wings during the 2010 outburst of V407 Cyg. 
These authors suggested that the emission and the asym- 
metric line profiles could originate in the wake produced by 
the convergent shock on the rear side of the companion star . 
We note that the scenario proposed by IShore et all (|201lh 
is compatible with the findings of the present paper. 
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